Character of electronic states in graphene antidot lattices: 
Flat bands and spatial localization 
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Graphene antidot lattices have recently been proposed as a new breed of graphene-based superlat- 
tice structures. We study electronic properties of triangular antidot lattices, with emphasis on the 
occurrence of dispersionless (flat) bands and the ensuing electron localization. Apart from strictly 
flat bands at zero energy (Fermi level), whose existence is closely related to the bipartite lattice 
structure, we also find quasi-flat bands at low energies. We predict the real-space electron density 
profiles due to these localized states for a number of representative antidot lattices. We point out 
that the studied low-energy, localized states compete with states induced by the superlattice-scale 
defects in this system which have been proposed as hosts for electron spin qubits. Furthermore, 
we suggest that local moments formed in these midgap zero-energy states may be at the origin 
of a surprising saturation of the electron dephasing length observed in recent weak localization 
measurements in graphene antidot lattices. 

PACS numbers: 73.23.-b, 73.63.-b, 73.21. Cd, 71.70.Di 



I. INTRODUCTION 

Investigations of the electronic properties of graphene 
constitute a relatively new and thriving sub-area of 
condensed-matter research.i After the first successful 
fabrication of graphene monolayers^ - by means of a me- 
chanical exfoliation of graphite - this two-dimensional 
semimetallic material has ignited tremendous experi- 
mental^i^i^ and theoretical^'^'^'i^iii'i^ii^'i^'i^ interests. 
From a fundamental standpoint, one of the main incen- 
tives for studying graphene stems from emergent analo- 
gies between the low-energy physics of the material and 
relativistic quantum mechanics ] Graphene is just 
as interesting from a practical point of view: owing to 
its exceptional properties - the extremely high mobility, 
chemical inertness, atomic thickness, and easy control of 
charge carriers by applied gate voltages - it holds a great 
promise for a carbon-based "post-silicon" microelectron- 
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Aside from simple graphene monolayers, patterning 
of monolayer films by nanolithography methods^^ - al- 
lowing feature sizes as small as tens of nanometers - 
has led to the demonstration of nanostructures such 
as Hall barsrii2^i2& quantum dots^^ nanoribbonsr2i2Si2£ 
and Aharonov-Bohm interferometers Another family 
of graphene-based structures has recently been proposed 
— triangular superlattices of holes (antidots) cut in a 
graphene sheet, known as antidot lattices i^!^^ Unlike 
pristine graphene which is semimetallic, graphene anti- 
dot lattices are semiconducting, with a direct band gap 
that depends on the antidot size. Square antidot lattices 
have been studied experimentally quite recently, corrob- 
orating the existence of a transport gapi^i^^ In addition, 
the weak localization correction to the conductance and 
a surprising saturation of the electron dephasing length 
at the superlattice scale have been observed in these ex- 



periments. 

In a recent work,'^^ antidot lattices have been proposed 
as a platform for quantum computation, with defects in 
this system envisioned as hosts for electron spin qubits. 
While the proposal of Ref. 32 focuses on localized states 
due to defects on the superlattice scale (such as missing 
antidots), their counterparts on the lattice scale (such as 
vacancies or adatoms) , essentially unavoidable along the 
antidot edges of experimental samples, are known to give 
rise to midgap (bound) states as welli^i^ One may thus 
expect that the two types of defects compete. In addi- 
tion, midgap states caused by lattice-scale defects might 
provide a plausible explanation of the maximal dephasing 
length observed in the experiment of Ref. [351: for suffi- 
ciently large charging energies, such midgap states can 
host local (spin) moments at the antidot edges that are 
known to be a very effective source of electron dephas- 
ingi^ This motivates us to systematically study midgap 
states in graphene antidot lattices, with a focus on their 
spatial profile and the corresponding dispersionless (flat) 
bands. 

On bipartite lattices (such as graphene) which have 
an excess of atoms on one of the two sublattices, zero- 
energy states are expected on very general groundsi^i^ 
As a consequence, midgap states may exist even in per- 
fectly symmetric and periodic antidot lattices. This was 
put forward by Shima and Aoki^^ in a general symmetry- 
based classification of superhoneycomb systems (i.e., tri- 
angular superlattices based on an underlying honeycomb 
lattice). Generically, however, midgap states will be in- 
troduced by irregularities in the shape of the antidot 
lattice on the atomic scale, which are unavoidable in 
present-day experiments. We analyze the resulting band 
structure and the spatial profile of the corresponding 
wave functions for a number of representative antidot lat- 
tices. Furthermore, we determine low-energy tunneling 
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FIG. 1: Graphene antidot lattice: (a) Lattice structure de- 
scribed by the basis vectors ai and a2, with magnitude 
|ai| — \a.2\ = La%/3 ; (b) hexagonal unit cell with a circular 
antidot; and (c) hexagonal unit cell with a triangular antidot. 
Vectors Si, S2, S3 specify positions of the nearest neighbors 
of a carbon atom on sublattice A. 



current distributions that can be compared with scan- 
ning tunnehng microscopy (STM) measurements. Our 
estimates for the charging energies of the predicted zero- 
energy states indeed suggest that for typical experimental 
parameters such states can form local moments and thus 
provide strong dephasing upon electron scattering from 
the antidot edges. 

The remainder of the paper is organized as follows. 
In Sec. nil we first briefly introduce the graphene su- 
perlattices of interest, accompanied by the notation 
and conventions to be used throughout (Sec. Ill Ap . and 
then lay out the framework for calculating their band 
structure (Sec. IIIB[) and the tunneling current distri- 
bution (Sec. IIICp . The obtained results are presented 
and discussed in view of the generic properties of the 
superhoneycomb- and bipartite systems in Sec. IIIII Fi- 
nally, we summarize our findings and conclude in Sec. lIVI 



II. DESCRIPTION OF ANTIDOT LATTICES 

A. Structure and nomenclature 

To set the stage, in this section we introduce graphene 
antidot lattices. A segment of a typical antidot lattice 
with a circular perforation is depicted in Fig. [TJa), with 
lattice basis vectors denoted by ai and a2. Its unit cell 
is a hexagon with an antidot in the center [Fig. IHb)]. 
We characterize the structure by the dimensionless side 
length of the hexagonal unit cell (L) and the radius of 
the antidot (R), both expressed in units of the graphene 
lattice constant a — 2.46 A. (Note that while L is an inte- 
ger, R can also take non-integer values; a = Occ'v/S, where 
ttcc = 1.42 A is the distance between nearest- neighbor 
carbon atoms.) Therefore, we use the notation {L, R} to 
specify antidot lattices with circular perforations. The 
unit cell of an antidot lattice with a triangular perfora- 
tion, which can be characterized by {L, D} (where D is 
the side length of the triangle), is shown in Fig.[T]^c). The 



number of carbon atoms (hereafter C atoms) per unit cell 
of an antidot lattice will henceforth be denoted by A^c- 

If we take a C atom on sublattice A at the ori- 
gin, its nearest neighbors are determined by vectors 
Si = (\/3/2, -1/2) flee, 62 = (0,1) flee, and <53 = 
(-V3/2,-l/2) acc (see Fig. p. Alternatively, for a C 
atom on sublattice B at origin, the corresponding vec- 
tors are —5i,—52, and — ^3. 



B. Method for band-structure calculation 

Given the large size of the unit cells in our supcrlat- 
tice - that in the cases of practical interest, with an- 
tidot diameter larger than about 10 nm, have numbers 
of C atoms in excess of several thousands - a band- 
structure calculation using the standard ab-initio meth- 
ods based on the density functional theory (DFT) is not 
feasible. We thus compute the electronic band struc- 
ture of graphene antidot lattices within a 7r-orbital tight- 
binding model. This method is known to reproduce very 
accurately the low-energy part of the DFT band struc- 
ture of pristine graphene 4^ 

The tight-binding Hamiltonian of an antidot lattice 
reads 
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(1) 



where vectors R designate unit cells (in total N of them) , 
d,„ (m = 1, . . . , Nq) specify positions of C atoms within 
a unit cell, S stands for nearest neighbors of a C atom at 
position R -I- dm, and t ~ 2.8 eV is the nearest-neighbor 
hopping matrix element, while the factor 1/2 is needed 
to correct for double-counting. After Fourier transforma- 
tion to momentum-space, the band-structure calculation 
amounts to a sparse-matrix diagonalization problem of 
dimension Nq. 

The Bloch wave functions corresponding to 
the energy eigenvalues SnaO^) are given by 
'0nka(r) = I]„ Cm''"Vmk(r), whcrc (jjmkir) = 
N^^/'^ X]r e*'^'^(/3(r — R — d„j). Here n enumerates en- 
ergy bands (with possible additional degeneracy labelled 
by a) and (^(r) denotes the 2pz orbital of a C atom. The 
energy eigenvalue problem i?|?/'k) = ^iV'k) reduces to 
det [H(k) — f S(k)] = 0, with matrices H and S given by 

Hmm'{^) = {4>mk\H\(j}m'k) and Smm'O^) = (</'mk 1 0m'k) , 

respectively. In the nearest-neighbor approximation 
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where the summation runs over superlattice vectors R = 
±ai, ±a2, ±(ai — a2) (only neighboring unit cells con- 
tribute), and (5(r,r') = 1 if C atoms at positions r and r' 



are nearest neighbors, while S 



(r,r'> 



otherwise. To a 



good approximation, the overlap of 2pz orbitals on differ- 
ent C atoms can be neglected, so that Smm'O^) = 5mm'- 
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This is a standard practice in the analyses of 7r-electron 
systems.'*^ 

For later reference, we describe a construction of an or- 
thonormal eigenbasis {V'nka} of H . To that end, we first 
note that 0„ik form an orthonormal set: ('/>mk|'/'m'k') = 
<5mm"5kk'- [For k = k' the orthogonality follows from 
S'mm'(k) = 5mm', whilc for k 7^ k' it holds be- 
cause 0mk and 0m' k' belong to different eigensubspaces 
of the lattice-translation operator.] This implies that 
(V'nkal'/'n'k'Q') — for k ^ k'. Therefore, to construct 
an orthonormal set, we find an orthonormal basis in the 
degenerate eigensubspaces of H(k). The orthogonality of 

the latter basis '^^{Cm^"''^)*Cm^'^ ^ — 5aa' implies that 

(V'nkalV'nk'a') = (^kk'l^aa'- 

C. Tunneling current distribution 

STM and related techniques are known as powerful 
tools for mapping out the spatial form of surface elec- 
tron states.^* In order to visualize the spatial structure 
of the studied low-energy states and provide the means of 
comparing our results to STM measurements, we predict 
the tunneling current distribution. 

The tunneling current 

rSp+eV 

J(r)cx / d£p{v-£), (3) 

J£f 

with £p being the Fermi energy and V the STM-tip bias 
voltage, can be used to probe the spatial dependence of 
the local density of states^ 

p(r;£) = ^|V'„ka(r)P<5[£-£:„(k)]. (4) 

nka 

This is a basis-independent quantity, expressed here in 
terms of an orthonormal eigenbasis {V'nka} of H, con- 
structed as described above. In the case that only one 
flat band at £ = falls into the energy window in 
Eq. ([3]), we have 

I(r)«El^"oka(r)P 

= E E Ck(r)<^™'k(r) . (5) 

ka mm' 

Strictly speaking, to calculate the tunneling current 
we would need to use the explicit form of 2pz orbitals. 
However, Eq. ([5]) simplifies if we assume that the 2^^ 
orbitals are well-localized on C atoms and hence neglect 
the overlap of the neighboring orbitals. After coarse- 
graining /(r) = (1/Ay) c^'^r' I(r') in the vicinity 
of a C atom positioned at r = R -I- d™, we obtain the 
sought-after tunneling current distribution 

/(r)«5]|C(r''''")p, (6) 

ka 

which is a lattice-periodic quantity. 




FIG. 2: Band structure for a {9,3} antidot lattice. Only 
bands above the Fermi level {£ = 0) are shown because of the 
particle-hole symmetry. F, K, M stand for the high-symmetry 
points in the Brillouin zone. 



III. RESULTS AND DISCUSSION 

In this section, we investigate the band structure of 
antidot lattices, with emphasis on localized states at low 
energy. The spatial profile of these states is characterized 
by a tunneling current distribution. We present results 
for both ideal lattices and lattices with defects (vacancies 
and/or debonded C atoms). Since we describe the system 
by a nearest-neighbor tight-binding model on a bipartite 
lattice, the resulting energy spectrum has particle-hole 
symmetry.— With this in mind, in the following we focus 
on the bands above or at the Fermi level (£ = 0). 

We start with an antidot lattice that is made up of 
perfect circular perforations, such as those considered in 
Ref. "35. The band-structure of an {L = 9, i? = 3} an- 
tidot lattice is shown in Fig. O A band gap f g opens 
at the F-point (k = 0). For antidot lattices with a rel- 
atively small number of removed C atoms (A'icm) com- 
pared to the total number of atoms in the original unit 
cell (iVtotai = ^rom + Nc), the band gap has been pre- 
dicted to scale as^ £g cx VA^rcm/A^totai- It has been 
demonstrated that localized states with energy inside the 
gap can be induced in this system by defects on the su- 
perlattice scale, e.g., missing circular antidots in the lat- 
ticei^ Such localized states have been proposed as hosts 
for local spin moments that may be utilized for quan- 
tum computation. In the following, we show that these 
midgap states appear generically in graphene antidot lat- 
tices, even without superlattice-scale defects. 

In the nearest-neighbor approximation, graphene has 
a bipartite lattice, that is, a lattice that can be divided 
into two sublattices, A and B, where only sites at differ- 
ent sublattices are connected through nonzero hopping 
matrix elements. Inui et al^ showed that in such sys- 
tems one has Na — Nb zero-energy states, where Na, Nb 
are the total number of sites on the respective sublat- 
tices. [This result was, in fact, implicitly known even 
earlier: it was obtained by Lieb"^^ as a prerequisite for 
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FIG. 3: Band structure of an antidot lattice with triangular 
antidots. Only bands above the Fermi level {£ = 0) are shown 
because of the particle-hole symmetry. F, K, M stand for the 
high-symmetry points in the Brillouin zone. The flat band 
at £ = is (a) sixfold degenerate (D = 6), and (b) ninefold 
degenerate {D = 9). 



the proof that the total spin in the exact ferromagnetic 
ground state of the Hubbard model on a bipartite lattice 
is S' = {Na — Nb)/2.] In a graphene superlattice with ua 
and ub sites per unit cell on the A- and B-sublattices, 
respectively, these states form ua — nn dispersionless 
(flat) bands at £ = Q. A sublattice imbalance ua — ?^B 
is generically introduced at the edges of graphene-based 
structures. Perhaps the most well-known example of the 
corresponding zero energy states are the edge states in 
zigzag graphene ribbons that form a partially flat band 
at the Dirac point ^li^ A pair of essentially flat (spin- 
polarized) bands close to the Fermi level was found also 
in a hydrogenated graphene ribbon, as demonstrated by 
ab-initio electronic structure calculations.**^ 

In antidot lattices one expects flat bands at zero energy 
due to sublattice imbalances along the edges of the anti- 
dots. The sublattice imbalance can occur even for perfect 
regularly shaped antidots. As an example, we study an- 
tidot lattices with triangular perforations (cf. Fig. [T]), 
which invariably have an imbalance of ua — ub — D 
per antidot. Consequently, a D-fold degenerate flat band 
emerges at 5 = 0, as depicted in Fig. [31 On the other 
hand, it is easy to check that antidot lattices with per- 
fect circular perforations always have ua = ub- These 
lattices therefore do not exhibit flat bands at f = (cf. 
Fig.HD. 

As pointed out by Inui et al.^ the single-particle 
states corresponding to the zero energy flat bands are 
pseudospin-polarized - they occupy only sites belonging 
to a particular sublattice. Figure 2] shows for an antidot 
lattice {L — lO^D — 9} how this effect manifests itself 
in the tunneling current distribution, which is propor- 
tional to the on-site electron density. Closer inspection 
reveals that the single-particle states ■(/'k corresponding 
to the flat band at f = indeed leave the B-sublattice 
sites unoccupied (amplitudes Cb^ = 0), while the A- 
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FIG. 4: Illustration of the spatial localization of the single- 
particle wave function corresponding to a flat band, (a) Tun- 
neling current distribution /(r), reflecting the population of 
only the A-sublattice sites. The sizes of the circles and the 
color brightness are proportional to /(r), which is scaled to 
< 7 < 1. The electron amplitudes obey the sum rule shown 
in the inset, (b) Tunneling current maxima along the y di- 
rection, at positions indicated by the arrows in (a). Inset: 
tunneling current distribution in the x direction. 



sublattice sites are occupied with amplitudes CAm (nor- 
malized as IC'a^P = 1), adding up to zero around 
each B-sublattice site [see Fig. HJa)]. Moreover, the elec- 
tronic states corresponding to the zero-energy flat bands 
are predominantly localized in the vicinity of the antidot 
edgci^ The extent of their spatial localization is illus- 
trated in Fig.IHb). 

It is important to emphasize that the above results are 
valid for an arbitrary system of bipartite structure, either 
infinite or of finite size, and even in the presence of off- 
diagonal disorder and/or an external magnetic field. 
(In the presence of a magnetic field, the hopping ma- 
trix elements become complex through the conventional 
Peierls substitution^* and the particle-hole symmetry is 
destroyed.) Therefore, even without actual calculation, 
we can conclude that the discussed flat bands at £" = in 
graphene antidot lattices remain flat in a magnetic field. 
Analogous results for finite-size graphene antidot flakes 
have recently been obtained numericallyi^ 
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FIG. 5: Effect of an on-site impurity potential Ve = —0.15 t 
along the antidot edges: (a) band structure, and the tunneling 
current distributions corresponding to (b) the lowest-lying flat 
band and (c) the highest flat band. 



FIG. 6; Effects of debonded C atoms in the antidot lattice 
{9,3.2}: (a) band structure, and the tunnehng current dis- 
tributions corresponding to (b) the low-lying quasi-flat bands 
and (c) the energy bands within the range 0.125 t < £ < 
0.257 t. 



A different perspective on the problem is furnished 
by a symmetry-based classification of superhoneycomb 
systems, as put forward by Shima and Aoki.'*^ They 
showed that such systems can have either semiconducting 
(with direct band gap), semimetallic, or metallic charac- 
ter. This classification turns out to depend not only on 
the global superlattice symmetry, but also on the specific 
atomic configuration within the unit cell. In particular, 
structures with 6m, 6m -I- 2, 6m + 3, and 6m -f 5 (m is an 
integer) atoms per unit cell belong to symmetry classes 
designated by Aq, Ac, Bq, and Be, respectively^i The 
respective degeneracies of the flat bands at f = in these 
classes are 61, 6/±2, 61 — 3, and 6/±l (/ > is an integer). 
According to this classification, our antidot lattices with 
circular perforations belong to the Ao-type of superhon- 
eycomb systems with I — 0, which is consistent with the 
absence of flat bands at £ = 0. 

In realistic graphene flakes, the dangling bonds along 
the edges are hydrogen-passivated, giving rise to an on- 
site potential. An on-site potential can also appear in a 
pure carbon system, due to a weak edge-magnetization 
induced by electron-electron interactions.^ Since the elec- 
tronic states corresponding to the zero-energy flat bands 
are largely localized at the antidot edges, one expects 
them to be strongly affected by such an edge potential. 
This motivates us to study the influence of an edge poten- 
tial, hereafter denoted by Ve, on the zero-energy states. 
We assume that Ve takes nonzero values at C atoms that 
have only two nearest neighbors. Such a potential breaks 
the particle-hole symmetry of the energy spectrum by 



destroying the bipartite lattice structure, since it couples 
sites that belong to the same sublattice. However, the 
main effect of Ve is a lifting of the degeneracy of the zero- 
energy flat bands, while their dispersionless character is 
largely being preservedi^ At the same time, the effect of 
this potential on the other bands is relatively small [cf. 
Figs. [5^ a) andlS^b)]. This finding is indeed analogous 
to the behavior of flat bands in nanoribbons with hy- 
drogenated edges, obtained using realistic first-principles 
DFT calculations.^- We stress that the observed effect of 
the on-site potential is robust, i.e., changing the magni- 
tude of Ve does not alter our qualitative conclusions. 

In what follows, we point out another generic feature 
of graphene antidot lattices: the occurrence of essentially 
flat bands at nonzero energies (even in the absence of 
an edge potential, Ve = 0). A characteristic example is 
presented in Fig. [51 While flat bands at f = in bi- 
partite lattices arise due to a global sublattice imbalance 
(riA ^ ub), as discussed above, we ascribe the quasi-flat 
ones at £■ ^ to local sublattice imbalances (while glob- 
ally UA = tib)- Such local imbalances can be induced 
even in regularly shaped antidot lattices, for instance, by 
debonded C atoms with a single neighbor at the edges 
[cf. Fig.llKb)]. 

The occurrence of quasi-flat bands may be understood 
as follows. A single debonded C atom induces a sublat- 
tice imbalance that leads to a zero energy "defect level" . 
One may view a collection of debonded C atoms as a col- 
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FIG. 8: Tunneling current distributions for a lattice with 
irregularly-shaped antidots, shown for the four lowest electron 
bands: f o = (brown), £i(k) (blue), f2(k) (red), and f3(k) 
(violet). The tunneling current magnitudes are represented 
by circles of different sizes. The sublattice polarization due 
to edge defects, along with the accompanying charge-density 
reconstruction, is shown in the zoomed part of the figure. 



FIG. 7: Effects of disorder in antidot lattices: (a) band struc- 
ture of an antidot lattice with defects on sublattice A (indi- 
cated by the arrow in the inset), (b) band structure in the 
presence of defects on both sublattices, and (c) unit cell with 
arrows indicating vacancies or C-adatom defects. The dashed 
curves in (a) represent the band structure of an ideal lattice 
(without disorder). F, K, M stand for the high-symmetry 
points in the Brillouin zone of the antidot lattice with the 
fourfold-enlarged unit cell. 



lection of "local sublattice imbalances" that induce one 
defect level each with wave functions that are localized 
in the vicinity of the defects. This picture is supported 
by the tunneling density of states shown in Fig.in^b). Lo- 
calized states induced by defects that are well separated 
from one another hybridize only weakly. Accordingly, 
these defect levels give rise to essentially dispersionless 
bands close to the Fermi level (£ = 0) As the distance 
between the defects decreases, the hybridization of de- 
fect states becomes stronger and the resulting bands are 
shifted away from the Fermi level. 

In disordered antidot lattices, debonded C atoms 
generically appear along the edges of irregularly shaped 
antidots, leading to a local sublattice imbalance. Such 
defects at the edge of a zigzag nanoribbon have been 
experimentally observed quite recently. In this experi- 
ment, the multilayer graphene structures were thermally- 
treated leading to edges that are mostly closed (i.e., 
folded from one layer to another). Open-edge structures 
with debonded C atoms are found in a local area where 
the folding edge is partially broken. Another source 



of local sublattice imbalances in disordered lattices are 
vacancies, for example, due to removed C atoms or C 
atoms that are sp^ rehybridized as a result of hydrogen 
chemisorptioni^i^i^ In Fig.[71[a), we show that a single 
defect in an otherwise perfectly circular antidot lattice 
produces, as expected, a zero-energy band. Needless to 
say, in realistic disordered systems defects do not repeat 
with the lattice period as assumed in Fig. [7(a). Some 
qualitative features of realistic disordered lattices can be 
captured by considering a superlattice with an increased 
unit cell. In Fig. [7]^b), we show a band structure for 
an antidot lattice with a unit cell containing four anti- 
dots with defects at different locations; this band struc- 
ture indicates that well-separated defect states weakly 
hybridize, giving rise to two pairs of quasi-flat bands that 
are symmetric with respect to the Fermi level. This sug- 
gests that for realistic disordered lattices one has a quasi- 
continuum of such low-energy states due to hybridization 
of defect states. 

We finally consider an antidot lattice with strong shape 
irregularities, as is likely the case in current experi- 
ments, '^^''^^ Although our calculation, the result of which 
is shown in Fig. [8l is based on periodically repeated de- 
fects, we expect that it will give reliable estimates for the 
localization properties of the appearing midgap states. 
The localization length found in this manner can be used 
to estimate the charging energy in the system. The ob- 
tained band structure has the expected features: the 
number of flat bands at f = is in agreement with gen- 
eral result for the bipartite systems, and the quasi-flat 
bands at nonzero energies correspond to local sublattice 
imbalances. The electron states for both types of flat 
bands are localized in the vicinity of the zigzag-like seg- 
ments of the antidot edges and debonded C atoms, as 
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illustrated in Fig. [8] (In contrast to the zigzag-like seg- 
ments, their armchair-like counterparts on the antidot 
edges do not exhibit electron localization.) Moreover, 
the tunneling current is maximal at the midpoints of the 
zigzag-like segments on the antidot edges. These features 
are in agreement with the existing experimental results 
on stepped graphite surfaces, which have edges with 
similar irregularities. 

Also depicted in Fig. [8] (see the zoomed part) is 
a "charge-density reconstruction" from the original 
graphene honeycomb lattice into a lattice with \/3 times 
larger period, which is rotated by an angle of 30°. This is 
a manifestation of a Friedel-oscillation-like phenomenon, 
that is, a long-range electronic perturbation caused by 
the presence of defects. Such interference phenomena, 
which are essentially a consequence of the wave nature 
of electrons, are familiar in the context of impurities on 
surfaces of metals®^ and were observed on the graphite 
and graphene surfaces using STM.— In the case of 
graphene, the charge-density reconstruction results from 
intervalley coupling of the electronic 7r-states. We thus 
conclude that there is strong intervalley scattering at the 
antidot edges in the lattice shown in F'ig. El This is in 
agreement with the experiments of Refs.[3Jand[3^, where 
the observed weak localization (instead of weak antilo- 
calization) is a signature of intervalley scattering at the 
antidot edges. 

The numerically obtained localization of the induced 
midgap states at the antidot edges, with a localization 
length on the order of a few interatomic distances, sug- 
gests that the charging energies for these states can 
be substantial. Assuming completely random edges of 
length Lg, one expects \nA — ~ \/ L^/a low-energy 
states per unit cell due to local sublattice imbalances 
(in a realistic disordered system these states will not be 
strictly at zero energy since the local sublattice imbal- 
ances UA — TLB in individual unit cells will partially can- 
cel between unit cells). For the parameters in the ex- 
periment of Ref. [ssl . one thus expects a spacing between 
these localized states of about \JLf.a ~ 10 nm along the 
edge, implying a charging energy that is substantially 
larger than the band gap in the system, which is roughly 
inversely proportional to the distance between antidots. 
Also, since the low-energy edge states due to local sub- 
lattice imbalances are separated by the antidot distance, 
one expects that their energy after hybridization is at 
most on the order of the band gap. Figure [7] supports 
this estimate. We thus conclude that the localized states 
observed in Fig. [8] have kinetic energies inside the band 
gap of the perfectly regular lattice, but charging energies 
that by far exceed that band gap. This suggests that 
local spin moments may form at the edges of disordered 
antidot lattices at Fermi energies where charge transport 



takes place. Since local magnetic moments are known to 
be an effective source of electron dephasing in weak local- 
ization experiments,— the studied midgap states offer an 
alternative explanation for the saturation of the electron 
dephasing length at a scale corresponding to the distance 
between antidots reported in Ref. l35|. 



IV. CONCLUSIONS 

In this work, we have studied the salient features of the 
low-energy band structure of graphene antidot lattices. 
Apart from strictly zero-energy (midgap) flat bands that 
arise from the bipartite lattice structure and a global sub- 
lattice imbalance, we have also found quasi-flat bands at 
low, but nonzero, energies that can be ascribed to local 
sublattice imbalances. In addition, we have examined the 
influence of an edge potential on the flat bands, showing 
that such a potential lifts the degeneracy of these bands, 
without affecting significantly their dispersionless char- 
acter. 

We have also investigated the spatial profile of the elec- 
tronic states corresponding to both classes of low-energy 
bands (flat and quasi-flat). By analyzing the tunneling- 
current distributions that can be compared to STM mea- 
surements, we have demonstrated that these electronic 
states are generically localized at the antidot-edges. The 
computed tunneling current distributions also show a 
charge-density reconstruction from the original honey- 
comb lattice to a lattice with -s/S times larger period and 
rotated through an angle of 30°. This phenomenon is 
indicative of intervalley scattering off irregular antidot 
edges, as also observed in recent experimentsi^!^ 

The spatial profiles of the localized midgap states that 
we find allow for a rough estimate of their charging en- 
ergies. That estimate suggests that such states can host 
local magnetic moments. We propose that such magnetic 
moments may be at the origin of a recently observed sat- 
uration of the dephasing length in graphene antidot lat- 
tices. "^^ In addition, the investigated midgap states com- 
pete with the localized states due to superlattice-scale 
defects, and therefore can have significant implications 
for a recent proposal of spin qubits in graphene antidot 
lattices.— 
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